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Abstract 

The isospin dependence of the spin-orbit interaction in hght neu- 
tron rich nuclei is investigated in the framework of relativistic mean 
field theory. The magnitude of the spin-orbit potential is consid- 
erably reduced in drip line nuclei, resulting in smaller energy split- 
tings between spin-orbit partners. The effect does not depend on the 
parametrization of the effective Lagrangian. The results are compared 
with corresponding calculations in the non-relativistic Skyrme model. 

The spin-orbit interaction plays a central role in the physics of nuclear 
structure. It is rooted in the basis of the nuclear shell model, where its 
inclusion is essential in order to reproduce the experimentally established 
magic numbers. In non-relativistic models based on the mean field approx- 
imation, the spin-orbit potential is included in a phenomenological way. Of 
course such an ansatz introduces an additional parameter, the strength of 
the spin-orbit interaction. The value of this parameter is usually adjusted 
to the experimental spin-orbit splittings in spherical nuclei, for example ^^O. 
On the other hand, in the relativistic framework the nucleons are described 
as Dirac spinors. This means that in the relativistic description of the nu- 
clear many-body problem, the spin-orbit interaction arises naturally from the 
Dirac-Lorenz structure of the effective Lagrangian. No additional strength 
parameter is necessary, and relativistic models reproduce the empirical spin- 
orbit splittings. 

Many properties of nuclei along the line of beta stability have been suc- 
cessfully described in the framework of models based on the mean-field ap- 
proximation. Conventional non-relativistic models that include density de- 
pendent interactions with finite range (Gogny) or zero-range (Skyrme) 
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forces have been extensively used to describe the structure of stable nu- 
clei. More recently, it has been shown that models based on the relativistic 
mean-field theory ^ provide an elegant and economical framework, in 
which properties of nuclear matter and finite nuclei, as well as the dynam- 
ics of heavy- ion collisions, can be calculated (for a recent review see [^). In 
comparison with conventional non-relativistic approaches, relativistic models 
explicitly include mesonic degrees of freedom and describe the nucleons as 
Dirac particles. Non-relativistic models and the relativistic mean-field theory 
predict very similar results for many properties of beta stable nuclei. How- 
ever, cases have been found where the non-relativistic description of nuclear 
structure fails. An example is the anomalous kink in the isotope shifts of Pb 
nuclei 0. This phenomenon could not be explained neither by the Skyrme 
model, nor by the Gogny approach. Nevertheless, it is reproduced very nat- 
urally in relativistic mean-field calculations. A more careful analysis has 
shown that the origin of this discrepancy is the isospin dependence of the 
spin-orbit term. With a spin-orbit term modified in such a way that it is 
similar to that derived in the relativistic mean-field model, the Skyrme model 
produces comparable results for the isotope shifts 0. Recently, in Ref 
another modification in the spin-orbit term of the energy functional of the 
Skyrme model has been proposed. An additional parameter was introduced, 
which is adjusted to reproduce the kink in the isotope shifts. Compared to 
conventional Skyrme forces, this approach produces a very different isospin 
dependence of the spin-orbit potential. 

Experiments with radioctive nuclear beams provide the opportunity to 
study nuclei with large neutron excess. Neutron drip lines of relatively light 
nuclei have become accesible, and the investigation of properties of such ex- 
otic objects is becoming one of the most exciting challenges in the physics 
of nuclear structure. Because of their relevance to the r-process in nucle- 
osynthesis, nuclei close to the neutron drip line are also very important in 
nuclear astrophysics. Knowledge of their stucture and properties would help 
the determination of astrophysical conditions for the formation of neutron 
rich stable isotopes |TD[. For drip line isotopes, nuclear shell effects become 



very important and the spin-orbit term plays an essential role. Very different 
scenarios for the isospin dependence of the spin-orbit interaction have been 
suggested in the Skyrme model |TT[, and in the relativistic mean-field the- 



ory |T^. In the present work we investigate the behaviour of the spin-orbit 
potential in light neutron rich nuclei. The description of drip line nuclei is 
complicated by the closeness of the Fermi level to continuum states. Pairing 
correlations, and the coupling between bound states and positive energy par- 
ticle continuum, are described in the Relativistic Hartree-Bogoliubov (RHB) 
model in coordinate space. 

In the Hartree approximation for the self-consistent mean field, the RHB 
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equations read 



A* -hj, + m + Xj WkJ ^"[vj- 

where ho is the single- nucleoli Dirac hamiltonian, and m is the nucleon 
mass. Uk and Vk are quasi-particle Dirac spinors, and Ek denotes the quasi- 
particle energies. The RHB equations are non-linear integro-differential equa- 
tions. They have to be solved self-consistently, with potentials determined 
in the mean-field approximation from solutions of Klein- Gordon equations 
for mesons. In the particle-particle {pp) channel the pairing interaction is 
approximated by a two-body finite range Gogny interaction 

VPP{1,2) = ^e-^'i-'^/^')' (W^i + B.P'' - H,P^ - MiP^P^), (2) 

1,2 

with parameters /ij, Wj, 5,, Hi and Mj {i = 1,2). 

In the present investigation we consider the isotopic chains of the even- 
even Ne and Mg nuclei. Systematic RHB calculations of ground state prop- 
erties have been performed, and the results reveal many interesting features. 
In this report we focus on the behaviour of the spin-orbit term. In the rel- 
ativistic mean-field approximation, the spin-orbit potential originates from 
the addition of two large fields: the field of the vector mesons (short range 
repulsion), and the scalar field of the sigma meson (intermediate attraction). 
In the first order approximation, and assuming spherical symmetry, the spin 
orbit term can be written as 

Vs.o. = --^VUr), (3) 
r or 



where Vis is the spin-orbit potential P, 113 



TT) 

Vis = (V-S). (4) 

^eff 

V and S denote the repulsive vector and the attractive scalar potentials, 
respectively, rrie/f is the effective mass 

meff = m-^{V - S). (5) 

In the following we present results for the even-even Ne and Mg isotopes. 
For the mean-field Lagrangian the NL3 [14] parametrization has been used. 



and the parameter set DIS [T^ for the finite range pairing interaction. Using 
the vector and scalar potentials from the self-consistent ground-state solu- 
tions, we have computed from (|^) - (|^) the corresponding spin-orbit terms 
for Ne and Mg chains. They are displayed in the upper panels of Fig. 1, as 
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function of the radial distance from the center of the nucleus. The magni- 
tude of the spin-orbit term Vs.o. decreases as we add more neutrons, i.e. more 
units of isospin. The reduction for nuclei close to the neutron drip is ~ 40% 
in the surface region, as compared to values which corespond to beta stable 
nuclei. This implies a significant weakening of the spin-orbit interaction. The 
minimum of Vs ^. is also shifted outwards, and this indicates the large spatial 
extension of the scalar and vector densities, which become very diffused on 
the surface. The main contribution to the densities in this region comes from 
the outermost neutron orbitals, whose wave functions are extened in space. 
The same effect can be observed in Fig. 2, where the rms radii are plotted 
as function of neutron number for Ne (a), and Mg (b) isotopes. We display 
neutron and proton rms radii, and the N^/^ curves. The dashed curves are 
normalized so that they coincide with neutron rms radii for ^ONe and ^'^Mg, 
respectively. Close to the neutron drip a sharp increase of neutron radii is 
observed, as compared to the mean-field curves 

This sudden mcrease 
indicates a possible formation of multi-neutron halo |]TB . 



In Fig. 3 we display the spin-orbit splittings of the neutron levels 

Els = En^lj=l-i/2 — -£'n,/,j=/+l/2, (6) 

for Ne and Mg isotopes as function of the neutron number. The neutron spin- 
orbit splittings decrease with neutron number. This is consistent with the 
gradual weakening of the spin-orbit term that is shown in Fig. 1. However, 
the result is at variance with non-relativistic mean-field studies that have 
used Skyrme forces [^. In the Skyrme model, the spin-orbit term included 
in the self-consistent mean field is of the form 

Vl' = W^{r){pxa), (7) 

with 

W.(r) = W.Vpr + W2Vpr,^r (8) 

where pr is the density for neutrons or protons (r = n or p) and r' is the 
opposite isospin. Wi and W2 are parameters. In the non-relativistic reduction 
of the Dirac equation the resulting spin orbit term takes the same form. 
However, the corresponding quantities Wi and W2 are coordinate dependent, 
and also introduce an isospin dependence which has its origin in the coupling 
constant Qp of the p- meson 



WiW = ^^SCl + Cl + Cl) (9) 



with Cf = {mgi/miY for i = a, a;, p. Non linear self-interaction terms for the 
sigma meson are included in the derivation of Wi and W2- 
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With the self-consistent proton and neutron densities of the drip hne 
nucleus ^°Ne, we have computed the neutron spin-orbit terms both in the 
Skyrme model and in the relativistic mean-field model. Spherical symmetry 
is assumed, and for simplicity, pairing has been neglected in this illustrative 
calculations. The effective force Skp has been used in Skyrme calcula- 
tions, and the NL3 parameter set for the effective relativistic Lagrangian. 
The results are summarized in Fig. 4. In Figs. 4a,b,d, and e, the neutron 
and proton densities, and the corresponding gradients are shown. The spin- 
orbit term calculated with the Skyrme force Skp is displayed in Fig. 4g (solid 
line), and the one resulting from relativistic calculations in Fig. 4i. The min- 
imum in the surface region is for the Skyrme force Skp «i 40% deeper than 
in the relativistic calculations. In order to have a more complete compari- 
son, in Fig. 4f we also display results of the relativistic model with the NLl 
0, [T|] (dotted line), and NL-SH (dashed line) effective forces. They 
turn out to be quite similar to those obtained with NL3. The small radial 
shift between the two minima is caused by the different isospin properties 
of the two parameter sets. The key point, however, is that in all relativistic 
calculations the spin-orbit term is much weaker compared to results of the 
Skyrme model. This of course explains the difference in the calculated spin- 
obit splittings. The spin-orbit term in the Skyrme model does not decrease 
in magnitude even for extremely neutron rich nuclei. This results in large 
spin-obit sphttings and quenching of the shell effects. 

The Skyrme model uses in principle an isospin independent two-body 
spin-orbit force. The exchange term to this force, however, induces a strong 
isospin dependence for the spin-orbit term in the self-consistent mean field. 
The resulting ratio for the W, parameters is: Wi/W2=2. On the other hand, 
in the relativistic approach the spin-orbit term is to a large extent a pure 
single-particle effect, as can be shown by a non- relativistic reduction of the 
Dirac equation. This is also the case in relativistic Hartree-Fock calculations. 
Therefore in the relativistic description there is no contribution from a two- 
body spin-orbit interaction, and consequently no exchange term. The isospin 
dependence of the parameters Wi and W2 in the relativistic theory comes 
from the p-meson. The contribution of the p mean-field increases with neu- 
tron number. This effect is illustrated in the lower panel in Fig. 1, where we 
plot the self-consistent p-meson potential for Ne and Mg isotopes. However, 
its contribution is much smaller than that of the a and uj fields, and there- 
fore the isospin dependence of Wi and W2 should be rather weak. In Figure 
4c we display the relativistic Wi and W2 calculated with the effective force 
NL3. The radial dependence comes from the effective mass m*. The two 
curves are very similar and in the surface region, where the spin-orbit term 
is peaked, their ratio is approximately W1/W2 ~ 1.1. If we now use the same 
ratio in the Skyrme model, the resulting spin-orbit term for the Skp force is 
displayed in Fig. 4g (dotted line). A significant reduction of the spin-orbit 
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term is observed, and the results are comparable to those obtained in the 
relativistic calculations. We have also computed the spin-orbit term for the 
effective interaction SkM* [21| and the two values of the ratio W1/W2. The 
results displayed in Fig. 4h show the same reduction in the surface region. 

In conclusion, we have shown that, in the framework of relativistic mean 
field theory, the magnitude of the spin-orbit potential is considerably re- 
duced in light drip line nuclei. With the increase of the neutron number, 
the effective one-body spin-orbit interaction becomes weaker. This result 
in a reduction of the energy splittings between spin-orbit partners. The re- 
duction of the spin-orbit potential is especially pronounced in the surface 
region, and does not depend on a particular parameter set used for the ef- 
fective Lagrangian. These results are at variance with those calculated with 
the non-relativistic Skyrme model. It has been shown that the differences 
have their origin in the isospin dependence of the spin-orbit terms in the two 
models. If the spin-orbit term of the Skyrme model is modified in such a 
way that it does not depend so strongly on the isospin, the reduction of the 
spin-orbit potential is comparable to that observed in relativistic mean-field 
calculations. 

This work has been supported in part by the Bundesministerium fiir 
Bildung und Forschung under contract 06 TM 875. 
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Figure Captions 

• Fig. 1 Radial dependence of the spin-orbit potential in self-consistent 
solutions for the ground-states of Ne (a), and Mg (b) isotopes. In the 
lower panel the self-consistent p-meson potentials are displayed. The 
NL3 parametrization has been used for the mean-field Lagrangian, and 
the parameter set DIS for the pairing interaction. 

• Fig. 2 Calculated rms radii for Ne (a), and Mg (b) isotopes as func- 
tions of neutron number. 

• Fig. 3 Energy splittings between spin-orbit partners, for neutron levels 
in Ne and Mg nuclei, as functions of neutron number. 

• Fig. 4 Self-consistent proton and neutron densities, the corresponding 
gradients, and the spin-orbit terms for ^°Ne. Results are displayed for 
the relativistic mean-field model and non-relativistic Skyrme calcula- 
tions. For description, see text. 
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Fig. 2 



